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Improving the Hadronization of QCD currents in TAUOLA and 
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We present our study of the hadronization structure of both vector and axial-vector currents leading to decays 
of the tau into two kaons and a pion. The cornerstones of our framework are the large- Nc limit of QCD, the 
chiral structure exhibited at low energies and the proper asymptotic behaviour, ruled by QCD, that is demanded 
to the associated form factors. The couplings of the theory are mostly constrained by this procedure and upon 
the analysis of BaBar data on e + e~~ — » KK-k we are able to predict the hadronic spectra. 



1. INTRODUCTION 

We focus on the nonperturbative dynamics of the 
strong interaction that can be learned thanks to 
the clean separation of the electroweak and QCD 
parts in the hadronic decays of the r lepton. Ad- 
vances in the study of the two meson processes 
have been presented in this workshop [I], [2]. 
The decay amplitude for the considered decays 
may be written as: 



M = --4= 



Gf 

V2 ud/us 



u Vr 7^(1 - ^)u T n^ 



(1) 



where our lack of knowledge of the precise 
hadronization mechanism is encoded in the 
hadronic vector, H u : 



(2) 



Symmetries help us to decompose Ti^ depend- 
ing on the number of final-state pseudoscalar (P) 
mesons, n. For three mesons in the final state, 
this reads: 

= V 1 ^(Q 2 ,s 1 , s 2 ) + V 2 ^(Q 2 ,s u s 2 ) + 
Q^(Q 2 , Sl ,s 2 ) + i V 3ll F^{Q 2 , Sl ,s 2 ) , (3) 



*I thank Stephan Narison and his collaborators for the 
organization of the QCD08 Conference. This work has 
been supported in part by the EU MRTN-CT-2006-035482 
(FLAVIAnet), by MEC (Spain) under grant FPA2007- 
60323 and by the Spanish Consolider-Ingenio 2010 Pro- 
gram CPAN (CSD2007-00042). 



and 



V x » = 

V 2 , = 

v 3 „ = 

Q» = 



(P2~Pl) V , 
(P3-PlY , 



£^uga P\ P 2 P3 ; 



(pi + P2 + Pz)n , Si = (Q- p l f . (4) 



Fi, i = 1,2, 3, correspond to the axial-vector cur- 
rent (-4^) while F4 drives the vector current (V^). 
The form factors F\ and F 2 have a transverse 
structure in the total hadron momenta, Q M , and 
drive a J p = 1 + transition. The scalar form fac- 
tor, Fs, vanishes with the mass of the Goldstone 
bosons (chiral limit) and, accordingly, gives a tiny 
contribution. This is as far as we can go without 
model assumptions, that is, we are still not able 
to derive the Fi from QCD. Our claim is that the 
most adequate way out is the use of a phenomeno- 
logically motivated theory that, in the energy re- 
gion spanned by tau decays, resembles QCD as 
much as possible. 

In any case, the approximate symmetries of QCD 
are useful. They rule what is the theory to be 
used in its very low-energy domain and guide the 
construction of higher-energy theories. 
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2. STRONG INTERACTION AMONG 
LIGHT-FLAVOURED MESONS 

Weinberg's Theorem [3] leads to conclude that 
a description in terms of the relevant degrees of 
freedom in the range of energies spanned by tau 
decays that incorporates the (approximate) sym- 
metries of QCD therein will be valid. Some phe- 
nomenological wisdom would tell that this proce- 
dure will be most adequate, indeed. Well below 
the p(770) mass, %PT [3], [4] is the effective the- 
ory of QCD. Being M T ~ 1.8 GeV, one needs to 
include the resonances as explicit degrees of free- 
dom. According to vector meson dominance [5], 
those of spin one would prevail. This is a result 
derived within Resonance Chiral Theory (RxT) 
[BJ worked in the large-A c limit of QCD. RxT 
is built upon the approximate chiral symmetry 
of low-energy QCD for the lightest pseudoscalar 
mesons and unitary symmetry for the resonances. 
The perturbative expansion of RxT is guided by 
its large- Nc limit |7J. At LO in 1/Nc, meson 
dynamics is described by tree level diagrams ob- 
tained from an effective local Lagrangian includ- 
ing the interactions among an infinite number of 
stable resonances. We include finite resonance 
widths (a NLO effect in this expansion) within 
our framework [8]. We are also departing from 
the Nc — > oo limit because we consider just one 
multiplet of resonances per set of quantum num- 
bers and not the infinite tower predicted that can- 
not be included in a model independent way. For 
convenience, we choose to represent the spin-one 
mesons in the antisymmetric tensor formalism [9] . 
The relevant part of the RxT Lagrangian is [6J, 

M- M- M- 
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VPPP i 



(5) 



where all couplings are real, being F the pion de- 
cay constant in the chiral limit. The notation is 
that of Ref. [BJ. P stands for the lightest pseu- 
doscalar mesons and A and V for the (axial)- 
vector mesons. Furthermore, all couplings in the 
second line are defined to be dimensionless. The 
subindex of the operators stands for the kind of 
vertex described, i.e., O vppp gives a coupling be- 
tween one Vector and three Pseudoscalars. For 
the explicit form of the operators in the last line, 
see [Ha, p], [TJ. 

Symmetries of low-energy QCD are not enough 
for RxT to describe it up to the perturbative re- 
gion. It is mandatory to match order parameters 
of spontaneous chiral symmetry breaking with 
partonic QCD related quantities. The matching 
of 2 and 3-point Green Functions in the OPE 
of QCD and in RxT has been shown to be a 
fruitful procedure [2], [10], p], [13], [14], [15]. 
Additionally, we will demand to the vector and 
axial-vector form factors a Brodsky-Lepage-like 
behaviour 16J. 

While symmetry fully determines the structure of 
the operators, it is the QCD-ruled short-distance 
behaviour who restricts certain combinations of 
couplings rendering RxT predictive provided we 
fix these few remaining parameters restoring to 
phenomenology. The set of couplings determined 
by this procedure depends on the particular pro- 
cess under study. In particular some couplings 
that were fitted phenomenologically in p° — * 
e + e - ,7r + 7r~ and a\ — > 7T7,p7r are predicted now. 

3. ASYMPTOTIC BEHAVIOUR AND 
QCD CONSTRAINTS 

There are 24 unknown couplings in C r x t (0) that 
may appear in the calculation of three meson de- 
cays of the t. They are reduced when computing 
the Feynman diagrams involved. There only ap- 
pear Fv, Gy, three combinations of the {Ai}f =1 , 
four of the {ci}J =1 , two of the {di}f =1 and four of 
the {(7i}f =1 - The number of free parameters has 
been reduced from 24 to 15. 
We require the form factors of the A 11 and V M cur- 
rents into KKtt modes vanish at infinite transfer 
of momentum. As a result, we obtain constraints 
[T2] , [T7] among all axial- vector current couplings 



but Ao, that are also the most general ones sat- 
isfying the demanded asymptotic behaviour in 
t — * 3irv T , studied along these lines in [15] . Pro- 
ceeding analogously with the vector current form 
factor results in five additional restrictions [12] . 
From the 24 initially free couplings in Eq. ([5]), 
only five remain free: C4, c\ + C2 + 8C3 — C5, 
di + 8 e?2 — d% , g± and 55. After fitting T(u> — > 3ir) 
-using some of the relations in Ref. only C4 

and 34 remain unknown. We obtained them from 
BaBar data on e + e _ — > KKir. 

4. RESULTS AND CONCLUSIONS 

Relating the measured isovector component of 
e + e~ — > KKtt to the total isovector cross-section, 
as done in Ref. [19] . is not trivial at all [20 . 
The issue is discussed in detail in our article [12] . 
We conclude that the measured isovector com- 
ponent in e + e~ — > KgK^n^ does not provide 
the total isovector component for the (partially) 
inclusive decay e + e~ — > KKtt. Once some as- 
sumptions are used, one can employ CVC to re- 
late e + e~ — > KKtt to r — > KKttv t . We fit- 
ted our expressions to BaBar data [21] obtaining 
c 4 = -0.044 ± 0.005. Our results for the decay 
widths of the considered channels are consistent 
with the PDG values [22 . To confront our pre- 
dictions for the spectra with the forthcoming ex- 
perimental data is a necessary task. 
We have obtained several sets of allowed values 
for those quantities that were still free after the 
high-energy conditions were imposed [12] . We get 
two predictions: 

• The ratio of the vector current to all con- 
tributions is, for all charge modes: 



T v (r -> KKtxv t ) 
r (r KKtxv t ) 



0.5. 



(6) 



• The spectra of the two independent KKtt 
decay modes of the r. We show in Fig- 
ure Q] that one for r — > K + K~ti~v t . The 
shape of the spectral function is similar for 
t — * K~ K®tt°v t . In both cases, the axial- 
vector current dominates the low-energy re- 
gion, while the vector one peaks at higher 
values of Q 2 . 



The plot presented here corresponds to: 

M ai = 1.17 GeV, c 4 = -0.04, g A = -0.5. (7) 



X --> K K It V 




Figure 1. Spectral function of r — > 
tained with the set of parameters (|7J. 



K- 



ob- 



We conclude by emphasizing that the phe- 
nomenological approach that we follow includes 
every relevant (and known) piece of QCD: We 
preserve the correct chiral limit at low energies 
that gives the right normalization to our form fac- 
tors, we employ large- Nc QCD arguments to use 
our theory in terms of mesons and it fulfills the 
high-energy conditions of the fundamental theory 
at the mesonic level. While we remain predictive 
by including only the lightest multiplet of axial- 
vector and vector resonances, our Lagrangian can 
be enlarged in a systematic way to account for the 
contribution of higher resonance states, whose pa- 
rameters may be fitted to experiment. 
This framework can easily be extended to study 
exclusive hadronic channels in e + e~ collisions, as 
we plan to do. Once the spectra for r — > KKttu t 
are presented, our work may shed some light on 
the issue of the importance of isospin violating 
corrections when relating tau and e + e~ hadronic 



cross sections. 

Our study can be very useful for the B-factories 
BaBar and Belle, that are already analyzing 
huge samples with good quality data, the tau- 
charm factory BES-III and a future super-B 
factory; specially, if our expressions are imple- 
mented in the currently used libraries for both 
types of processes, namely TAUOLA [23] and 
PHOKHARA [21]. 
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